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Abstract 


this work we calculate the mass spectrum, weak decay constants, two photon decay widths, and two 

I 

P-igluon decay widths of ground (IS), and radially excited (2S, 3S,...) states of pseudoscalar charmoniuum 
£jand bottomonium such as pc and rjb, as well as the mass spectrum and leptonic decay constants of ground 
CO state (IS), excited (2S, ID, 3S, 2D, 4S,...,5D) states of vector charmonium and bottomonium such as J/b, 
and T, using the formulation of Bethe-Salpeter equation under covariant Instantaneous Ansatz (CIA). Our 
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O mass spectra as well as the transition amplitudes. However, the price we have to pay is to solve a coupled 
77 set of equations for both pseudoscalar and vector quarkonia, which we have explicitly shown get decoupled 


results are in good agreement with data (where ever available) and other models. In this framework, from 
the beginning, we employ a 4 x 4 representation for two-body {qq) BS amplitude for calculating both the 


in the heavy-quark approximation, leading to mass spectral equation with analytical solutions for both 
H masses, as well as eigenfunctions for all the above states, in an approximate harmonic oscillator basis. 
The analytical forms of eigen functions for ground and excited states so obtained are used to evaluate the 
decay constants and decay widths for different processes. 


X 


1 Introduction 

Heavy flavor mesons are of the type QQ, or Qq, in which one of the quarks belongs to a heavy flavour 
such as c or b. Quarkonia usually refers to heavy mesons with quark composition {QQ). Charmonium and 
bottomonium are two types of quarkonia that are observed with quark composition cc and bb respectively. 
Physically, the charmonium system is analogous to positronium system. However, historically, its role as 
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a model system in QCD is similar in importance to the role of Hydrogen atom in Quantum mechanics, 
due to the fact that this system is non-relativistic and has a hydrogen atom like spectrum. Charmonium 
in fact occupies a valuable intermediate position within QCD, being neither in the purely non-relativistic 
regime nor the regime where chiral symmetry breaking dominates the physics. This makes it a relatively 
clean system in which to study non-perturbative QCD dynamics, QCD-inspired quark-potential models as 
well as lattice QCD, which have been rather successful in describing the observed features of the spectrum 
(see [1] for a review). However, charmonium cannot be considered to be completely understood; as an 
example, in recent years a number of new charmonium resonances have been claimed in experiment, several 
of which cannot be easily reconciled with the predictions of simple quark-potential models. Charmonium 
spectrum provided a simple example of how QCD works, which was made even more compelling with the 
subsequent observation of bottomonium spectra. Investigation of properties of these mesons gives a lot of 
information about heavy quark dynamics. Below the heavy flavour meson pair (DD) production threshold, 
all charmonia states have been well established, and there is a good agreement between predictions based 
on potential models and data. However, there are many charmonium and charmonium like states [21 Emils] 
observed above DD threshold in the past ten years. Some of these are good candidates for the charmonia 
predicted in different models. However, many have unusual quantum numbers [21 El IH E] (which can not be 
reached by pure quark states), indicating that exotic states such as multi - quark states, molecule, hybrid, 
or glueball may have been observed. The spectroscopy and decay rates of quarkonia are quite important 
to study as huge amount of high precession data acquired from many experimental facilities world over 
are continuously providing accurate information about hadrons particularly in charm and beauty flavour 
sectors. The mass spectrum and the decays of these bound states can be tested experimentally and 
thus studies on them may throw valuable insight of the heavy quark dynamics and lead to a deeper 
understanding of QCD further. 

As regards the dynamical framework, to investigate these properties is concerned, many non-perturbative 
approaches, such as Lattice QCD Chiral perturbation theory [9], heavy quark effective theory 

[To], QCD sum rules mn N.R.QCD [13], dynamical-equation based approaches like Schwinger-Dyson 
equation and Bethe-Salpeter equation (BSE) [TTl [151 EHl 113 EEl EHl 120], and potential models [211 122] 
have been proposed to deal with the long distance property of QCD . Some of the interesting works on 
heavy-quark spectroscopy in the framework of BSE are [231121]. 

Many theoretical predictions on the decay properties particularly the leptonic and two-gamma de¬ 
cays of quarkonia based on the relativistic quark model or potential model [2S112SII2Z], Bethe-Salpeter 
equation |281129] . heavy-quark spin symmetry PI and lattice QCD ED are available in literature. Such 
studies have become a hot topic in recent years, due to observation of many new states at various high 
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energy accelerators at BABAR, Belle, CLEO and BES-III collaborations [32l [33l [SU |35]. All this has 
opened up new challenges in theoretical understanding of heavy hadrons and provide an important tool 
for exploring the structure of these simplest bound states in QCD and for studying the non-perturbative 
(long distance) behavior of strong interactions. 

We wish to mention that though heavy quarkonia (cc and bb) are very well described in terms of 
NR potential models, and relativized potential models (like [22]) can treat light mesons simultaneously 
with heavy mesons, however, they cannot be regarded as an ideal theory that is closest to QCD, due 
to their inability in implementing Lorentz and gauge invariance for a consistent treatment of motion 
of quarks within the hadron. Most of such models, employ the corrections, or the replacement 

—)■ (p^ + in the kinetic energy operator (for each quark) in a Schrodinger type equation, which 

provides the bulk of momentum dependent effects. However, interaction terms require a more substantial 
form of treatment consistent with relativity. We further wish to point out that though, in heavy quarkonia, 
NR treatment may at best be a good hrst approximation, but in charmonium (where charmed quark is not 
heavy enough to be considered non-relativistically), one hnds typical velocities, )( ~ 0.4, making relativistic 
effects in charmonium important, which is specially so for electroweak decay properties as well as various 
other transitions at high energies, as mentioned above. These decay widths are sensitive to relativistic 
effects (see [IS])- Thus NR and semi-relativistic equations (as in [22]) are not completely adequate as 
proper starting point for light, and semi-light quark systems. 

In this respect, the Bethe-Salpeter Equation (BSE), with its 4D Lorentz-invariant structure is a natural 
candidate for hadron physics, since it is hrmly rooted in held theory, and its general framework rests on 
Lorentz and gauge invariance, due to which it has a very wide applicability, all the way from low energy 
spectroscopy to high energy transition amplitudes, and can be used to study both light and heavy quark 
systems in an integrated manner. Thus the BSE framework provides a realistic description for analyzing 
hadrons as composite objects. Despite its drawback of having to input mo del-dependent kernel, these 
studies have become an interesting topic in recent years, since calculations have shown that BSE frame¬ 
work using phenomenological potentials can give satisfactory results on more and more data are being 
accumulated. The Lorentz-invariance at the input form of BSE formally ensures a continued validity of 
predictive powers from low energy spectroscopy to high energy processes within a common dynamical 
framework. The 4D BS wave function [JH] IM] ET] |3El EH] HO] (that carries the entire non-perturbative 
information) arises as solution of the 4D BS equation. Since this wave function is determined as a solution 
of the formal dynamical equation, it does not suffer from the uncertainties of a variational determination 
(as in NR approaches [22]), and is able to carry over the micro causality information contained in the BSE. 
From non-perturbative 4D BS wave function, one can work out the 4D hadron-quark vertices, that are 
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effective coupling vertices of hadron with all its constituents (quarks), and then naturally work out tran¬ 
sition amplitudes for different processes over a vast range of energies, employing Feynman diagrammatic 
techniques over various quark-loop diagrams corresponding to the process considered. 

A number of recent studies in the framework of BSE have been carried out on leptonic decays and 
two photon decays of mesons. The decay into two photons is considered as an interesting experimental 
playground in the mesonic physics of the near future. The two-photon decay of mesons can be used to 
identify the flavor of quark and anti-quark states. Two photon couplings also provide a useful probe of 
the internal structure of mesons. The states such as, rjc and rji, have been observed to have partial widths 
consistent with quark model predictions. These states can also decay into two gluons [H], which accounts 
for a substantial portion of the hadronic decays for states below the cc or bb threshold. It is clearly 
important to have accurate quark model predictions of widths for all experimentally accessible gg-mesons. 
Thus, in this work we are predominantly interested in studying the mass spectrum and leptonic decays of 
ground and excited states of pseudoscalar and vector quarkonia, such as J/-^, and T respectively 

(which proceed through the coupling of quark- anti-quark loop to the axial vector and vector currents 
respectively), as well as, the two-photon radiative decays, and two-gluon decays of ground state (IS") and 
radial excitations (25* and 3S) of conventional pseudoscalar mesons such as pc and rji, in the framework 
of BSE under Covariant Instantaneous Ansatz (CIA) [TSl 1201 EZl EHl EH SO]) though some of the excited 
states have not yet been experimentally discovered. 

We have employed the BSE framework under Covariant Instantaneous Ansatz (CIA). CIA is a Lorentz- 
invariant generalization of the instantaneous ansatz (lA). What distinguishes CIA from other three- 
dimensional (3D) reductions of BSE is its capacity for a two-way interconnection [TSl 1201 EH SO]: an 
exact 3D BSE reduction for a gg system (for calculation of the mass spectrum), and an equally exact 
reconstruction of original 4D BSE (for calculation of transition amplitudes as 4D quark loop integrals). In 
some of our previous works in BSE under CIA, we had made use of a conhnement potential V (g, (f) given 
in Eq.(33), for a number of successful predictions for not only the mass spectrum |TH ESI SD] of light (gg) 
and some heavy mesons, but also a number of processes such as leptonic decays of vector mesons such as: 
p, u, 0, J/'0, and T [10] through the process D —>■ e+ - 1 - e~, leptonic decays of pseudoscalar mesons such as 
TT, K, D, Ds and B through the process, P —)■ I -|-lg and two-photon decays of tt, and rjc mesons (22] through 
process, P —?• 7 -|- 7 , as well as the processes proceeding through quark triangle loop diagrams involving 
two or more hadron-quark vertices, such as: single photon decays of and uj^ mesons through process 
V —)■ P -|- 7 , as well as strong decays of hrst radial excitations of lightest vector meson states p(1450), 
and a;(1420)[12] through process V' ^ V + P (where V, and P refer to vector and pseudoscalar mesons). 
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However this form of confining potential was employed mostly for light mesons, and for some cases to 
heavy mesons. However, in our previous BSE framework [38l |39l HQ], we had adopted a 16 x 1 column 
representation for two-body {qq or qq) BS amplitude (though both 16 x 1 and 4x4 representations of 
BSE are completely equivalent Cl)- We had employed Gordon reduction on the mass shells of individual 
quarks to obtain the 3D mass spectral equation. And for bringing out the structure of the full 4D BS wave 
function in a 4 x 4 matrix form (that is needed for calculation of transition amplitudes through quark-loop 
diagrams), we made use of standard transformations (Hj of charge conjugation of spinors. This is in 
contrast to the BSE framework employed in this paper, where we have employed a 4 x 4 framework of 
BSE from the very beginning (with the same conhning potential in Eq.(33)), and used it for both the 
mass spectral predictions as well as the transition amplitude calculations involving heavy quarkonia (both 
pseudoscalar and vector), though the price one has to pay for this approach, is to solve a set of coupled 
equations, to obtain the masss spectral equations. However, in the present work, we have explicitly shown 
that in the heavy quark approximation (valid for cc, and bb systems), these equations can be decoupled, 
and analytical solutions (both mass spectrum and eigen functions) of these equations can be obtained 
using approximate harmonic oscillator basis. 
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Table 1: Quantum numbers of lowest quarkonium states. These quantum numbers are repeated for the 
radial excitations of these states 

Towards this end, we hrst analyze the quarkonium states according to their total angular momentum, 
Tf = ^ + parity P = (—l)'^"''^, and charge conjugation, C = (—1)'^+'^ for states classihed as with 
the lowest possible states for quarkonia are listed in Table 1. We start with the most general structure 
of Bethe-Salpeter wave function for = O”’*' and = 1 . It can be seen from Table 1, that 1 

state has not only the ground state (IS*) components, but also the orbital (ID) excitations. The same 
holds true for their radial excitations. We have started with the full Dirac structure of these states with 
all the Dirac covariants multiplying various scalar functions of internal hadron momentum, q as in H, 
with various Dirac structures incorporated into the wave functions in accordance with the power counting 
rule [Ml So] suggested recently. We then put the formulated BS wave functions into the instantaneous 
BSE and turn the equation into a set of proper coupled equations [13] for the components which appear 
in the formulation. These equations are then explicitly shown to decouple in the heavy-quark limit, and 
are reduced to a single mass spectral equation, whose analytic solutions in an approximate harmonic 
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oscillator basis yield the mass spectrum and wave functions for ground and excited states of rjc, rjb, J/ip 
and T. For this we need 6 input parameters (that include two input quarks masses me, and rrib). We then 
derive the leptonic decay constants, and the decay widths for two-photon decays, and two-gluon decays 
of pseudoscalar quarkonia rje, and rjb for their ground and radially excited states, as well as the leptonic 
decay constants for ground, radially and orbitally excited states of vector quarkonia, J/ip, and T, with 
the 6 input parameters hxed above. 

This paper is structured as follows: in section 2 we introduce the BSE and formulate the instantaneous 
BSE for pseudoscalar and vector quarkonia. In section 3 we start with the generalized formulation of BS 
wave functions for pseudoscalar = O”’*') and vector = 1 ) quarkonia, with dehnite quantum 

numbers and derive their mass spectral conpled eqnations. In section 4, we study the leptonic decays 
of pseudoscalar and vector qnarkonia. In section 5, we derive the two photon and two gluon decays of 
psendoscalar qnarkonia. Finally, section 6 is relegated to nnmerical resnlts and discussion. 


2 Formulation of BSE under CIA 


Lets consider a qq comprising of fermionic qnarks of masses mi and m 2 respectively. We start with a 4D 
BSE for qq system, written in a 4 x 4 representation of 4D BS wave fnnetion 4/(P, q) as: 

S?‘(Pi)4<(P,9)S?‘(-P2) = d'q'K[<,,q'mP,q') (1) 

where K{q,q') is the interaction kernel between the qnark and anti-quark, and pi ^2 are the momenta 
of the qnark and anti-qnark, which are related to the internal 4-momentnm q and total momentum P of 
hadron of mass M as. 


Pi,2fi = ± q^ 


( 2 ) 


/ 2 2 \ 

where mi ^2 = are the Wightman-Garding (WG) definitions of masses of individual quarks 

which ensure that on the mass shells {P.q = 0) of either qnarks, even when mi 7 ^ m 2 . However for eqnal 
mass mesons, (mi = m 2 = m), we have mi = m 2 = 


Then pi 2 ^ becomes, 

Pi,2p = ^Pf,±qbi 


(3) 


Now it is convenient to express the internal momentnm of the hadron q^ as the sum of two parts. 
They are: (i) the transverse component, q^ = q^i — {q ■ P)Pbi/P^ which is orthogonal to total hadron 
momentnm (ie. q ■ P = 0 regardless of whether the individual quarks are on-shell or off-shell), and 
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(ii) the longitudinal component, aP^ = {q ■ P)P^/P‘^, which is parallel to P^. Thus we can decompose 
as, = {Ma,q), where the transverse component, q is an effective 3D vector, while the longitudinal 
component, Ma plays the role of the time component. The 4-D volume element in this decomposition is, 
d'^q = d^qMda. To obtain the 3D BSE and the hadron-quark vertex, use an Ansatz on the BS kernel K 
in Eq. (1) which is assumed to depend on the 3D variables g^, g|, as, 

;r(g,g')=i^(g,g') (4) 

Hence, the longitudinal component, Ma of g^, does not appear in the form K{q,q') of the kernel. For 
reducing Eq.(l) to 3D form, we dehne 3D wave function 'ijj{q) as: 

m = ^J Mda^iP,q) (5) 

Substituting Eq.(5) in eq.(l), with definition of kernel in eq.(4), we get a covariant version of Salpeter 
equation, 

-mi)T(P,g)(^2+ "^ 2 ) = y ^^A:(g,g')^(g'), (6) 

and the 4D BS wave function becomes. 


= SF{pi)T{q)SF{-p2) 


(7) 


where 

( 8 ) 

plays the role of hadron-quark vertex function. The kernel is taken to have 7 ^ 07 ^ form, details of which are 
given in the next section. Thus with this form of the kernel, we can write r(g) = / (f„^D(g, g')7^V’(?07^ 
(where V is spatial part of the kernel). We can to a good approximation express 7/i'0('?O7^ ~ 
where 0 involves the spin-spin interactions alone, that factor out of the RHS of the hadron-quark vertex 
on taking the dominant Dirac structures in 'ip{q^) in the calculation of 0, and we can write r(g) = 
® ^ And, Sf{p) is the usual fermionic propagator of the quarks, given as. 


Sf{±Pi,2) 


^1,2 =^"^72 

Ap 2 


(9) 


where Ai ^2 = P 12 P 2 ) which can also be decomposed as in[T9] , 


with 


SF{±Pi) = 


A+(g1 


+ 


A* (g) 


I{i)Ma + — Ui I{i)Ma + ^M + uj^ 


( 10 ) 


f 

Ai^(g) = ^[^Wi±/(i)(mi + |)] (11) 
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where i = 1,2 for quark and anti-quark respectively, and I{i) = (—Here Af(g), are called as 
projection operators. With the projected wave functions, one can rewrite the BSE as, 

(1m t cui + Ma){^M T 0 J 2 - q) = A±(g)r(g)A±(g), (12) 

where projected wave functions, are obtained by the operation of projection operators on ' 0 (g) as, 

0^^(g) = A^(g)^0(g)^A^(g). 

With contour integration over da on both sides of Eq.(12), we can obtain: 

_ _AfWTOAj® K(q)T(mm 

M — (jJi — UJ2 M + UJi UJ2 

The complete wave function can separate in four parts as: 

0 (g) = 0 ++(g) + 0 +"(g) + 0 "+(g) + 0 ""(g). 

The BSE then reduces to four independent equations as [15] : 

(M - 2cn)0++(g) = -A+(g)r(g)A+(g) 

(M + 2 n;) 0 --(g) = Ah(g)r(g 0 A 2 -(g) 

0 +-(g) = 0 -+(g) = 0 

where ui = UJ 2 {= cn), for equal mass systems. In fact the four equations constitute an eigenvalue problem 
that should lead to evaluation of mass spectra of pseudoscalar (for preliminary work, see [13]) and vector 
charmonium and bottomonium states such as rji,, J/il: and T. The framework is quite general so far. 
Thus to obtain the mass spectral equation, we have to start with the above four equations to solve the 
instantaneous BS equation. 


(13) 

(14) 

(15) 

(16) 


3 Derivation of mass spectral equations of pseudoscalar and 
vector qucirkonia 

In this section, we show the method to solve an Instantaneous BS equation, which is hrst applied to the 
calculation of the mass spectrum of equal mass heavy pseudoscalar and vector quarkonia. We hrst write 
down the most general formulation of the relativistic BS wave functions, according to the total angular 
momentum (J), parity (P) and charge conjugation(C) of the concerned bound state. We then put this BS 
wave function in Eq.(16) to derive the mass spectral equations, which are a set of coupled equations. We 
now illustrate this procedure for pseudoscalar and vector quarkonia. 




• For pseudoscalar mesons, the complete decomposition of 4D BS wave function in terms of various 
Dirac structures and scalar functions (j)j{P,q) multiplying them is [14] : 


q) = {01 (P, q) + q) + <IUP. q) + [f. 0 ^{p. q)}i^ (i7) 

where (j)j = (f)j{q‘^,q.P,P^); (j = 1,2, 3,4) are the Lorentz scalar amplitudes multiplying the various 
Dirac structures in the BS wave function, \1'(P, g). Now in our framework, by use of a naive power 
counting rule in [3H1 EHl SO], we had shown that the Dirac structures associated with amplitudes 0 i 
and 02 are leading, while the structures associated with 03 and 04 are sub-leading, and would con¬ 
tribute lesser to meson observable calculations in comparison to leading Dirac structures associated 
with 01 and 02 - And in various calculations [HHl ESI HO], we had shown that among the two leading 
Dirac structures associated with amplitudes 0i, and 02 , the structure associated with 0i (i.e. 75 ) is 
dominant. We wish to point out that Munczek and Jain jH] have also earlier shown that 0i is the 
dominant amplitude for all ground state pseudoscalar mesons. V.Sauli |15] has recently shown that 
the same is valid for their excited states as well. And this is particularly true for mesons made from 
heavy flavour quarks. 

In the center of mass frame, where = (0,g), we can then write the general decomposition of the 
instantaneous BS wave function for pseudoscalar mesons {J^'" = 0 '"), of dimensionality M in the 
center of mass frame as [TH] : 

0^(g) [M0i(g) + /’02(g) + iMq) + ^Mq)h 5 , (is) 


where 0i, 02 , 03 and 04 are even functions of g and M is the mass of the bound state (of the 
corresponding meson). We now obtain the algebraic forms of these amplitudes. 


For this, we put the instantaneous BS wave function 0'^(g) (Eq.(18)) into the last two equations of 
Eq.(16), and this leads to independent constraints on the components for the Instantaneous BS wave 
function as in [IH1EII21]: 


= 0 ; 04 = 


-02M 


m 


(19) 


So we can apply the obtained constraints Eq.(19) to Eq.(18) and rewrite the relativistic wave function 
of state (O”’*') as. 


0-^(g) [M0i(g) + /’ 02 (g) + —02(g)]75- (20) 

m 

One can then see that the Instantaneous BS wave function of 0 state is determined by only two 
independent functions 0i, and 02. Putting the wave function Eq.(20) into the hrst two equations of 
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Eq.(16) and by evaluating trace over the y-matrices on both sides, we obtain two BS coupled integral 
equations, 


(M - 2u) 
(M + 2u) 


Mq ) + —(1 - 


m 




»i(« - -(1 - 

m m'* 


= epy 

= —0p 


(27r)3 
f d?q' 


( 27 r 


E(g,g' 
E(g,g') 


0i(g') + -(l + ^02(g') 


m 




—(1 + 


m 




To decouple these equations, we hrst add them. Then we subtract the second equation from the hrst 
equation. For a kernel that can be expressed as V{q — (f) = V{q)S^{q — ^), we get two algebraic 
equations which are still coupled. Then from one of the two equations so obtained, we eliminate 
01 (g) in terms of 02 (<?), and plug this expression for 0i(g) in the second equation of the coupled 
set so obtained to get a decoupled equation in 02(g)- Similarly, we eliminate 02(g) from the second 
equation of the set of coupled algebraic equations in terms of 0i(g), and plug it into the hrst equation 
to get a decoupled equation entirely in 0i(g). Thus, we get two identical decoupled equations, one 
entirely in 0i(g), and the other that is entirely in, 02(g)- Employing the limit, a; ~ m on RHS, these 
equations can be expressed as: 


M 2 

~T 

M 2 

~T 


02 ^ 

m- - q- 0i(g) = 0pE(g)m0i(g) + -^V (g)0i(g) 
-Tii^ -f 02(g) = QpV{q)m4>2{q) + ^E^(g)02(g)- 


( 22 ) 


It is to be mentioned that these two decoupled Eqs. (22) would resemble h.o. equations, but for 

V (g) term on the RHS of these equations. It will be shown later from the dehnition of kernel in 

_ 2 

Eq.(31), and Eq.(34), that this second term involving V (g) is negligible in comparison to the hrst 
term involving V (g) on the RHS, and can be dropped, and these equations would then resemble exact 
harmonic oscillator equations. However, due to identical nature of these equations, their solutions 
are written as: 0i(g) = 02 (g) ~ 4>p{q), which represent the eigenfunction of the pseudoscalar meson 
obtained by solving the full Salpeter equation and we can write the wave function for = 0 ^) 
state as: 

0^(g) ^ [M + + ^]0p(gl75- (23) 

m 

In view of the arguments made above, in the structure of 0'^(g) for pseudoscalar mesons, M 75 would 
have dominant contribution among all the other Dirac structures ( and this is more so for heavy QQ 
mesons) 

For vector mesons, the complete decomposition of 4D BS wave function in terms of various Dirac 
structures is 
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q) = ixi + ifX2 + [g ■ e - if\Xz + - <lf)]Xi 

(g-e)X5 + {q-e)^X6 + {q-e)^X7 + {q-P){f<^ - if)x^ (24) 

where Xa = Xa{q‘^,q-P,P‘^)]{c( = 1,2,3,..., 8 ) are the Lorentz scalar amplitudes multiplying the 
various Dirac structures in the BS wave function, \t'(P, q). We again mention that in our framework, 
by use of a naive power counting rule in [3811391110], we had shown that the Dirac structures associated 
with amplitudes Xi and X2 are of leading order. Those with amplitudes Xs, ...,X 6 are sub-leading. 
We had also shown that, the structures associated with X 7 and Xs are more much more suppressed 
than even the sub-leading Dirac structures, and would contribute very little to meson observable 
(specially heavy mesons) calculations in comparison to Dirac structures associated with xi,...,X 6 - 
Thus we ignore the Dirac structures associated with xt and Xs (see [101 for details). And in various 
calculations [Ml HQ], we had shown that among the two leading Dirac structures associated with 
amplitudes xi, and X 2 , the structure associated with xi he. y.e is most dominant. In Ref. mmsi, 
it was also shown that Xi is the dominant amplitude for not only ground state vector mesons, but 
also for their higher excitations. 

We thus write the instantaneous BS wave function ijj^{q) with dimensionality M for vector quarkonia, 
= 1 ) up to sub-leading order (i.e. 0{1/M^) as per the power counting scheme [10]) as in 

case of pseudoscalar mesons, in the center of mass frame as: 

?/>^(g) = M/xi + ifx2 + w - g-e]X3 + [fH - {.q-Pjf]^Xi + (g-e)X5 
+(9.«)/'§ (26) 

where xi, ••■,X 6 are even functions of q and M is the mass of the bound state (of the corresponding 
meson). We now derive the mass spectral coupled equations from Eq.(16). Putting Eq.(25) into the 
last two equations of Eq.(16), and we obtain the independent constraints on the components for the 
Instantaneous BS wave function: 

Mxi 

Xs =- 

m 

-X 2 M 

X4 = -, 

m 

(26) 

with Xe = Xs = 0- Applying the constraints in Eq.(26) to Eq.(26), we can rewrite the relativistic 
wave function of state (J^*" = 1 ) as: 

i^iq) = [Mi + q.i-]xM + Pt]x-2(q), (27) 

m mm 
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where, we have been able to express the instantaneous wave function in terms of only the 

leading Dirac structures associated with amplitudes Xi and X2- Putting the wave function Eq.(27) 
into the hrst two equations of Eq.(16) and by evaluating trace over the 7 -matrices on both sides, we 
can obtain two independent BS coupled integral equations: 

{M -2uj) [xi(g) - —X2{q)] = ©v / -l!^V{q,q) \xi{q') - —X2{q') 

L m \ J [iTiy L m \ 

{M + 2uj) [xi(g) + —X 2 (g)l = ©v / ^^V{q,q') [-Xi(g') “ —X 2 (g') 

L m \ J {2Tiy V m J 

(28) 

To decouple these equations, we proceed as in the above case of pseudoscalar mesons. We hrst add 
these above equations. Then subtract the second equation from the hrst. And for a kernel expressed 
as V{q — ^) = V{q)6^{q — ^), we get two algebraic coupled equations in xi and X 2 - We Eliminate xi 
from the hrst equation in terms of X2, and plug it in the second equation to get an equation entirely 
in terms of X2- Similarly, we eliminate X2 from the second equation in terms of Xi, and plugging in 
the hrst equation to get an equation entirely in terms of Xi- We thus get two identical decoupled 
equations, one in Xi(g), and the other in X 2 {q), and in the approximation, u ^ m, (valid for heavy 
quarkonia), these equations can be expressed as: 

^ -m^ -f xi(g) = ©vF(g)mxi(g) + 

X2{q) = QvV{q)mx2{q) + ^v‘^{q)X2{q), (29) 

It is to be mentioned that these two decoupled equations would again resemble h.o. equations, but 
_2 

for V (q), term on the RHS of these equations. It will be shown later from the dehnition of spatial 
part V of the kernel in Eq.(31), and Eq.(34), that this term is negligible in comparison to the hrst 
term involving V on the RHS, and can be dropped, and these equations would then resemble exact 
harmonic oscillator equations. 

From the identical nature of both these equations, we can again express the solutions as, Xi(q) = 
X 2 (q) ~ (t>v{.q)) and we can write the wave function for = 1 ) as: 

'ip^iq) ^ [M/ -h q.e— + ( 30 ) 

m mm 

In the structure of 'ip^{q) for vector meson, as argued above, is the most dominant Dirac structure, 
among all the other Dirac structures above (and this is more so for heavy QQ mesons). We wish to 
mention that the framework is quite general up to this point. We now introduce the BS kernel. 
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Now, as regards the BS kernel K{q, q') [15113811391110] is concerned, it is taken to be one-gluon-exchange 
like as regards the spin dependence ( 7 ^{S> 7 ^), and color dependence (^Ai.|A 2 ), and has a scalar part V: 


K{q,q') = q') 

3 f 1 

y{Q,Q ) = '^^qq J + 4mim2AoM>r^)“5 

wjg = 4:M^mim2ula,{Ml) 

fi,r2\ Ml . 

>1 

K = (1 -|- 4111 , 1111 , 2 AQM ‘^ r ‘^)~"2 

M> = Max{M, nil + 1 ^ 2 )■ 


a 


OJ, 




(31) 


The scalar part is purely conhning (as in Ref. [HI |23| H] , and the Martin potential [16] employed for heavy 
mesons). Here the proportionality of cjT on as{Q^) is needed to provide a more direct QCD motivation 
[T5] to conhnement, and Wg is postulated as a spring constant which is common to all flavors. Here in 
the expression for H(g, g'), the constant term Cq/uiI is designed to take account of the correct zero point 
energies, while the Hg term (Hg << 1 ) simulates an effect of an almost linear conhnement for heavy 
quark sectors (large mi, m 2 ), while retaining the harmonic form for light quark sectors (small mi, m 2 ), as 
is believed to be true for QCD (see HSlEHlSn]). Hence the term (1 + 4 mim 2 Ho^^^^) ^ in the above 
expression is responsible for effecting a smooth transition from harmonic {qq) to linear {QQ) conhnement. 


We now try to work on the spatial part V{q,(f) of the conhning potential K{q,(f). If we take the 
parameter Hg = 0 (which corresponds to case of light mesons {qq), since due to Hg << 0 , the square root 
factor in the denominator, k = (1 -f- 4 mim 2 HoA 4 ^^^"^)~^ = !)> and V would look like: 


Q 


V'{q,q') = y^, I d^r[r^-^] 


9_]pi(q-q')-f 


OJn 


(32) 


(where V' is V without the factor k in denominator). Making use of the fact that —^yis the fourier 


transform of in momentum space, and the colour factor for bound qq system, (|Ai. 4 A 2 ) = 


- 3 , we can 


write 


V(f,?) =75(27r)3|yj+ _5-), (33) 

Uq 

which can be written as, V'{q,(f) = V'{q)5^{q — q'). Putting this spatial part in place of the spatial 
part of the kernel V in Eq.(22) and Eq.(29) (with the second term on the RHS of each of these equations 
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dropped for reasons explained next), these eqnations wonld resemble harmonic oscillator eqnations, and 
we can have analytical solntions for them. Bnt for parameter Aq ^ 0 (i.e if we reintrodnce the factor k), 
we can write the complete potential V{q, (f) as: 

V{q,q') = V{q)5%q-q'), 

V{q) = + ^](27^)^ 

UIq 

K={1- (34) 

—2 — 

As mentioned above, we have dropped the terms with V , in comparison to the terms involving V on 

02 

the RHS of Eqs. (22) and (29), as the coefficients, ^py = -^^^qq associated with the former have a 
very small contribntion (< 0.638%) in comparison to the coefficient ffpy = mBpyCjA associated with the 
latter for both psendoscalar and vector mesons, dne to << cjA for rjc, rjb, J/il^ and T. The nnmerical 
valnes of these coefficients, and their percentage ratio for both psendoscalar (r/c, hfc)) vector (J/'0 and 
T) mesons are given in Table 2 below, which jnstifies these terms being dropped. 



Hp 

n'p 

fi/p% 

fly 

n'y 

nv 

Vc 

0.0558 

0.000356 

0.638% 




Vb 

0.4564 

0.00202 

0.4426% 




J/^ 




0.0246 

0.00006889 

0.280% 

T 




0.2293 

0.0005115 

0.2230% 


02 

Table 2: Nnmerical valnes of coefficients, Vtp = m0pa;A, and Vt'p = -^ujqq associated with the terms 
involving V and V respectively, for psendoscalar mesons rjc and rjb in RHS of Eqs. (22), and their percentage 
ratio, along with the corresponding valnes, and percentage ratio of Q,v = "^^v^qq, and ffy = for 

vector mesons, J/'tp and T on RHS of Eqs.(29). The inpnt parameters of our model are: Cq = 0.21, 
ojq = .ISGeR, QCD length scale A = 0.200GeR., Ag = 0.01, and the input quark masses, rric = 1.49GeR., 
and rrib = 5.070GeV. 

To derive the mass spectrum, we put the spatial part, V in Eq.(34) into the equations Eq.(22), and 
Eq.(29) for psendoscalar and vector mesons respectively. But the form of V in Eq.(34) suggests that 
these equations have to be solved numerically. However, to solve these equations with Ag % 0, we follow 
an analytical procedure on lines of Ref. [3H1I1Z1, where we treat k as a ’’correction” factor due to small 
value of parameter, Ag [Ag « 1), while we work in an approximate harmonic oscillator basis, due to 
its transparency in bringing out the dependence of mass spectral equations on the total quantum number 
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N. (In this connection, we wish to mention that recently, harmonic oscillator basis has also been widely 
employed to stndy heavy qnarkonia using a Light-front quark model [H].) The latter is achieved through 
the effective replacement, «:=(! — ^ (1 -|- 2Aq{N + |))~T (which is quite valid for heavy cc 

and bb systems), in Eq.(34) on lines of [36lll7j. With this, we can reduce Eq.(22) and Eq.(29) to equations 
of a simple quantum mechanical 3D- harmonic oscillator with coefficients depending on the hadron mass M 
and total quantum number N. The wave function satishes the 3D BSE for equal mass heavy pseudoscalar 
and vector mesons respectively as given below: 


and 









\J'\- + 2Aq[N + |) 


+ ^]4>p{q) 


(35) 



f)<t>v{q) = Qvmugg 




I -+ ^] 4 >v{q), 

^l + 2Ao{N+l) ^0 


(36) 


where, with use of dominant Dirac structures [39l SO], we can to a good approximation express, 0p = 
—4 for pseudoscalar states, and 0y = —2 for vector states, and write Eqs.(35) and (36) in the same 
expression as: 

Epy(j)py{q) = {-fipy^ q + q'^)4>py{q), (37) 

where, ~ ^ + f ^ ~ I = 0,1,2,... correspond to S,P,D,... states respectively. Epy = 

^ - Tip + + “^AoiN + |), and (j)py{q) are the eigen functions of Eq.(37) for equal mass heavy 

pseudoscalar, and vector qnarkonia. The inverse range parameters for pseudoscalar and vector meson 
respectively are, f3p = (4yj^=^^fy)T and (3v = and are dependent on the input 

kernel and contains the dynamical information, and which only differ from each other due to spin-spin 
interactions. We can express Eq.(37) for the I = 0, 2,..., states studied here as: 


(l>p,viq) + + -]^[Epy 

q Ppy 


+ l)Pf,y 


(p 


q^](t)py{P) — 0 . 


(38) 


q 

Assuming the form of the solutions of the above equation as (ppy{p = h{q)(fe 7^, where p = Ppy, 
we obtain the equation. 


,21 2 2g, ,E 21 3,,_, 

/! (9) + ( J J - ^)h (9) + 1^ - ^ - ^M‘l) = O' 


( 39 ) 


Making use of the power series method to solve this equation, we obtain the energy eigenvalues as: 
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( 40 ) 


{Epy)N — 2f3^{N + -); N — 2n + l;n — 0,1,2, 


with the normalized forms of energy eigen fnnctions for I = O(S'), and for I = 2{D) states derived as: 


0P{V')(15', q) — 


TT' 


1 

e 


3 / 4 ^ 3/2 


cf)py^{2S,q) - (2)^^^^3/4^3/2 

^vO-D,q) = (—) ! ^3/4^7/2 

(j)py){3S,q) = 


8 ^ 7r^/^/33/2 


_jL 

(l-3^)e 

g 2 g-^ 

(1 _ ^ + 


( l ) v { 2D , q ) (15^374^7/2(1 7^2)^^ 


15/?2 15^4 


4 

e 2/3% 


, r/i c _ . 35 , 4/2 1 ,4 _ ^ 84^4 

( t > P { v ){ , q ) (^g) ^3/4^372(1 105/32 + 

2 




(l>v{^D,q) ! ^3/4^7/2 




63/32 53^4 


105/34 105/36 

4 

)g 2 g 2/32 _ 


e 2.32 


(41) 


We will nse them for a description of the eqnal mass heavy psendoscalar and vector mesons. In Fig.l, 
and Fig.2, we now give the plots of these normalized wave fnnctions Vs. q (in Gev.) for different states 
of psendoscalar cc, and bb qnarkonia. And in Fig.3, and Fig.4, we give the plots of these wave fnnctions 
Vs. q (in Gev.) for different states of vector cc, and bb qnarkonia. It can be seen from these plots that the 
wave fnnctions corresponding to nS, and nD states have n — 1 nodes. 



Fignre 1: Plots of wave fnnctions for states (IS', ...jdS") Vs q (in Gev.) for psendoscalar cc states. 
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Figure 2: Plots of wave functions for states {IS, ...,45') Vs q (in Gev.) for pseudoscalar bb states. 


The mass spectrum of ground (IS) and excited states for equal mass heavy pseudoscalar (0” 
vector (1 ) mesons respectively is written as: 

1 M2 ' 


2 (— -™ + — 


and 


2/32 ^ 4 
1 M2 


^^\ll + 2Ao{N +l)) = N + %N = 2n +l;n = 0,1,2 ,..., 

u. 


2/32 ^ 4 


(^ - m2 + + 2Ao{N + ^))=N + ^-,N = 2n + l-,n = 0,l, 2.. 


, )2 


and 


(42) 


(43) 



BSE - GIA 

Expt.[35] 

Pot. Model pi] 

QGD sum rule [12] 

Lattice QGDjS] 

m 


2.9509 

2.983±0.0007 

2.980 

3.11T0.52 

3.292 

2.981 

^Vc(2S) 

3.7352 

3.639±0.0013 

3.600 


4.240 

3.635 

^»?c(3S) 

4.4486 


4.060 



3.989 

^»?c(4S) 

5.1048 


4.4554 



4.401 

^VbSS) 

9.0005 

9.398 ±0.0032 

9.390 

9.66± 1.65 

7.377 

9.398 

^Vb{2S) 

9.7215 

9.999±0.0028 

9.947 


8.202 

9.990 

MribiSS) 

10.4201 


10.291 



10.329 

Mr,b{iS) 

11.0968 





10.573 


Table 3: Masses of ground and radially excited states of rjc and rjh (in GeV.) in present calculation (BSE- 
GIA) along with experimental data, and their masses in other models. 
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BSE - CIA 

Expt.[35] 

Rel. Pot. Model [49] 

Pot. Model [2T] 

BSE[52] 

Lattice QCD[5T] 


3.0974 

3.0969± 0.000011 

3.096 

3.0969 

3.0969 

3.099 


3.6676 

3.6861± 0.00034 

3.685 

3.6890 

3.686 

3.653 


3.6676 

3.773± 0.00033 

3.783 


3.759 


Afp(35) 

4.1945 

4.03± 0.001 

4.039 

4.1407 

4.065 

4.099 


4.1945 

4.191±0.005 

4.150 


4.108 



4.6856 

4.421T0.004 

4.427 

4.5320 

4.344 



4.6856 


4.507 


4.371 


Afp(5S) 

5.1463 


4.837 

4.8841 

4.567 



5.1463 


4.857 




Mt(is) 

9.6719 

9.4603± 0.00026 

9.460 

9.4603 

9.460 


Mt(2S) 

10.1926 

10.0233±0.00031 

10.023 

9.9814 

10.029 


Mx(i_d) 

10.1926 


10.154 


10.139 


Mx(3S) 

10.6979 

10.3552±0.00005 

10.355 

10.3195 

10.379 


Mt(2D) 

10.6979 


10.435 


10.438 


Mt(4S) 

11.1887 

10.5794±0.0012 

10.586 

10.5995 

10.648 


Mr(3D) 

11.1887 


10.704 


10.690 


Mr{5S) 

11.6657 

10.865±0.008 

10.869 

10.8465 

10.868 


Mt(4D) 

11.6657 


10.949 




Mr{6S) 

12.1296 

11.019±0.008 

11.088 

11.0713 



AfT(5D) 

12.1296 







Table 4: Masses of ground, radially and orbitally excited states of heavy vector quarkonium, J/-^ and T 
in BSE-CIA along with their masses in other models and experimental data (all units are in GeV). 
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Figure 3: Plots of wave functions for states (15, 3D) Vs q (in Gev.) for vector cc states. 

The input parameters of our model are: Cq = 0.21, ujq = .ISGeV, QCD length scale A = 0.200GeV., 
Aq = 0.01, and the input quark masses, rric = 1.49GeV., and mb = 5.070GeV. The results of mass spectral 
predictions of heavy equal mass pseudoscalar and vector mesons for both ground and excited states with 
the above set of parameters is given in table 3 and 4. We now present the calculation of the leptonic decay 
constants fp and fy for these equal mass heavy pseudoscalar and vector mesons, rjc, r]b, J/4’ and T. 

4 Leptonic Decays of equal mass Heavy pseudoscalar and vec¬ 
tor quarkonia 

We now do the calculation of decay constants of equal mass pseudoscalar and vector mesons such as, r^c, 
rjb, J/^l^ and T, which are defined as [43l IHH] . 

ifpPfi =< 0|Q*7m75Q|-P > • 

fvMe.iP) 0\Q-f,Q\V{P) > (44) 

The decay constants fp and fy thus can be evaluated through the loop diagram which gives the coupling 
of two-quark loop to the axial vector current and vector current respectively, and can be expressed as 
a quark-loop integral (for some of our recent works on leptonic decays in the framework of BSE, see. 
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Figure 4: Plots of wave functions for states (15, 3D) Vs 



q (in Gev.) for vector bb states. 


[Ml Ea sni ss]): 

fpP^ = '/3j j^Trl>I''’(P,(/)np75|, 

fvM(, = 731 rr[4,>'(p, (45) 

where is polarization vector of vector meson satisfying e.P = 0. These equations can be reduced to 3D 
forms by dehning the 3D wave function, 'ip{q) = — f q). Thus we can write Eq.(45) as: 

fpPf, = Vs J ^^Tr[^/>-^(g)7^75] 

/yMe^ = Vs J ■^^Tr[ip^{q)'y^] (46) 

where (following Eq.(23) and Eq.(30)), the complete 3D Bethe-Salpeter wave function of state (O”’*') and 
(1 ) is rewritten as: 


tl^^iq) = Np(j)p{q)[M + ^ + —]75, 


m 

M 


i^^iq) = Nv(j)v{q)[M^ + q.e -h /Z’ + 

m m 


.e 


m 


(47) 


where Np and Ny are the standard BS normalizers which enter into the BS wave function, and m and 
M are the masses of quarks and the corresponding meson respectively. 

The 4D BS normalizer Np and Ny are evaluated from the current conservation condition: 

d 


2tP^ = 


(27) 


-Tr < 4/(P, q) 


dP. 


Sp {pi 


vl>(P,g)5^^(-p2) +(1-2). 


(48) 
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Carrying out derivatives of inverse quark propagators of constituent quarks with respect to total hadron 
momentum P^, evaluating trace over products of gamma matrices, following usual steps, we then express 
the above equation in terms of the integration variables q and a. Noting that the 4D volume element 
d^q = (PqMda, we then perform the contour integration in the complex a- plane by making use of the 
corresponding pole positions. Then integration over the variable q is hnally performed to extract out the 
numerical results for Np and Ny for different equal mass pseudoscalar and vector mesons. The above 
equation due to the orthogonality condition, P.q = 0 and P.e = 0 reduces to a simple forms given below 
for equal mass pseudoscalar and vector mesons: 

d^q (f 


Np^ = 


16M 


m 


and 


iVy^ = IQMm 


(27r)3 a; 
d?q 




4>v('i)- 


(49) 


(50) 


(27r)3 a;3 

Here, in the above equations for decay constants and the BS normalizers, the wave functions, (pp and 
(pv represent the eigenfunctions of the pseudoscalar and vector mesons obtained by solving the full Salpeter 
equation. Their algebraic expressions are taken from Eq. (41), and are then employed to calculate fp and 
fv as well as the BS normalizers Np, and Ny for ground as well as the excited states of rjcfih, J/tp, and 
T. The numerical values of BS normalizers of all these states are given in Table 5 below. 




N,. 

Nj/ip 

Np 

15 

5.2239 

5.7441 

5.9102 

6.4666 

25 

3.1632 

3.6679 

3.6788 

4.1743 

ID 



3.6233 

4.1534 

35 

2.3845 

2.8657 

2.8317 

3.2904 

2D 



2.8079 

3.2805 

45 

2.1531 

2.4103 

2.3591 

2.7881 

35) 



2.3462 

2.7058 


Table 5: Numerical values of BS normalizers for ground state and excited states of rjci Vb, J/'P T (in 
GeV units) in present calculation. 

The decay constant from Eq.(46) becomes: 

fpP„ = V3Np [ d^4>pmrl(M + t>+ ^)7j 
J (27r)3 m 

fvMe„ = VsNr f pP-4,v(q)Tr[(M I + q.e- + q> + P^- 

J (27r)3 m mm 

(51) 
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Evaluating trace over 7 -matrices, and carrying out integration over d^q, we obtain: 

fp = iVsNp J 

fv = 4:V3Nv J 



BSE - CIA 

Expt.fM] 

Lattice QCD [6] 

Pot. Model pr] 

QCD sum rule [12] 

fridlS) 

0.4044 

0.335±0.075 

0.3928 

0.471 

0.260T0.075 

fvci^S) 

0.3308 



0.374 


fvc{3S) 

0.2908 



0.332 


fvbiis) 

1.0168 


0.667 

0.834 

0.251T0.072 

frib{2S) 

0.8066 



0.567 


fvb(3S) 

0.7134 



0.508 



Table 6 : Leptonic decay constants, fp of ground state (IS) and excited state (2S) and (3S) of rjc and rjb (in 
GeV.) in present calculation (BSE-CIA) along with experimental data, and their masses in other models. 

The calculated values of decay constants are given in Tables 6 and 7. We next calculate the decay 
widths for two photon and two gluon decays of pseudoscalar quarkonia rjc and rjb and their radially excited 
states in the next section. 


5 Two photon and two gluon decays of pseudoscalar quarko- 
nium 




Figure 5: Diagrams contributing to process P — 77 . The second diagram is obtained from the hrst 
diagram by reversing the direction of internal fermion lines. 

The two-photon decays of these states have been the subject of numerous studies aimed at further 
understanding the accuracy of theoretical models of the charmonium and bottomonium systems based on 
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BSE-CIA Expt.[35] 

Pot. Model [21] 

BSE |52] 

LatticeQCD[50| 

Light front model [48] 

fj/ip(is) 

0.3745 

0.411T.007 

0.317 

0.459±.028 

0.399T 0.004 


fip{2S) 

0.2953 

0.279±.008 

0.253 

0.364±.024 

0.143± 0.081 

0.2474 

fipilD) 

0.2897 

0.210±0.00024 

0.243±.017 



Ui3S) 

0.2610 

0.174±.018 

0.226 

0.319T.022 



/p(2D) 

0.2348 

0.1424±0.0033 

0.157T.011 



/p(4S) 

0.2399 

0.1608±0.0016 

0.288±.018 




0.2570 

0.1424±0.0033 

0.157T.011 



fr(is) 

0.9005 

0.708±.008 

0.645 

0.498±.020 


0.1822 

fr{2S) 

0.6072 

0.482±.010 

0.439 

0.366T.027 


0.1944 

/t(1D) 

0.5972 



0.261T.021 



/t(35) 

0.5446 

0.346T.050 

0.393 

0.304±.027 



fr{2D) 

0.5755 



0.155±.011 



/t(45) 

0.5859 

0.3406T.00037 

0.259±.022 



/t(3D) 

0.6274 



0.155±.011 




Table 7: Leptonic decay constants, fy of ground state (IS) and excited state (2S) and (3S) heavy vector 
quarkonium, J/'ij) and T in BSE-CIA along with their masses in other models and experimental data (all 
units are in GeV). 
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the available data. The process, P —)■ 77 proceeds through the famous quark triangle diagrams shown in 
Fig.5. In this hgure the second diagram is obtained from the first one by reversing the directions of the 
internal fermion lines in the quark loop. Let, ki^ 2 i and ei ^2 be the momenta and polarization vectors of the 
two outgoing photons. Let, pi and p 2 be the momenta of constituent quark and anti-quark constituting 
the hadron with the total momentum P = pi + P 2 and relative momentum 'F(P, g) is the 4D hadron 
Bethe-Salpeter wave function. For sake of convenience, we introduce the relative momentum of the two 
outgoing photons: 2Q = ki — k 2 - In terms of P and Q, we can express the momenta of the outgoing 
photons as: ki ^2 = |P ± Q. The momenta of the third quark in the two diagrams can be expressed as: 
p = q — Q and p = q + Q respectively. 

The amplitude for the process in Fig.5 can be expressed as the sum of amplitudes for the two diagrams 
in this hgure as: 


My.(P^77) - Q)h + hSF{'J + Q)h]} 


(53) 


Cq = 36 for charm quark and Cq = -^e for bottom quark, Spiq T Q) are the quark and anti-quark 


= +2 

3 

propagators which are given as: 


SpiqPQ) = 


T0)±rn 


(54) 


{qTQyTm^' 

Now for heavy hadrons like cc and 66 , where the system is basically non-relativistic, it is convenient to take 
the internal momentum q « M, and hence, g^ << Q^, where it can be easily seen that Q'^ = ^- We 
now make use of the fact that the 4D volume element, d^q = dP’qMda, and use the relationship between 
3D and 4D BS wave functions. And since the rest of the integrand does not involve g, we can express the 
amplitude in the above Eq.(53) as: 


Mf,{P ^ 77 ) = - iQ)i,]}. (55) 

Now for a pseudoscalar meson with state = 0 the general relativistic Salpeter wave function 
'ip{q) can be written as: 

^^(g) = iVp0p(g)[M + f+ ^]75 (56) 

where Np is the standard BS normalizer which enters into the BS wave function, and m and M are the 

masses of quarks and the corresponding quarkonium respectively. The wave function 4>p{q), represents 

the eigenfunction of the Pseudoscalar meson obtained by solving the full instantaneous Bethe-Salpeter 

__F 

equation. For ground state (IS*) mesons, 0p(g) = e ^. The wave functions for excited states can 
similarly be expressed, while the BS normalizer Np is given as in Eq.(49). Pntting the Salpeter wave 
fnnction ip^iq) into the amplitnde in Eq.(55), and performing the trace over the gamma matrices, we 
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obtain: 


Mfi{P 77 ) = [Fp] 

^ fiua(jPfi^luQa^2a') (57) 

where Fp is the radiative decay constant for two-photon decays of psendoscalar meson and is expressed 


as: 


Fp = 




^ Np 


d?q 


M) 


m‘^ + M‘^/A J (27r)^ 

The decay width for P —)■ 77 is related to the corresponding decay constant, Fp by the expression: 

\Fp\^M^ 


r 


P^77 


bdvr 


(58) 


(59) 


The nnmerical resnlts of Tp_^^^ are given in the Table 8 . We now work ont the process of two-gluon decays 
of psendoscalar quarkonia. 



BSE-CIA 

Expt.pS] 

BSE|2n] 

BSEpE] 

RQM[2B] 

ra 

Pot.Model pT] 

rT;c(15)^77 

7.9178 

7.2±1.2 

7.14 

3.5 

5.5 

7.5-10 

11.17 

rr;c(25)->-77 

5.7889 


4.44 

1.38 

1.8 

3.5-4.5 

8.48 

rr;c(35)-5>77 

0.2995 



0.94 



7.57 

r)? 6 (lS)-s >77 

0.7376 


0.384 

0.22 

0.35 

0.560 

0.58 

^’? 6 ( 2 S )-!>77 

0.5076 


0.191 

0.11 

0.15 

0.269 

0.29 

i^ 7 ? 6 ( 3 S )->-77 

0.4261 



0.084 

0.10 

0.208 

0.24 


Table 8 : Two-photon decay widths of ground state (IS) and excited state (2S) and (3S) psendoscalar 
mesons, t]c and rji, in present calculation (BSE-CIA) along with their masses in other models and experi¬ 
mental data (all values are in units of Kev). 

The two-gluon decay width gives information on the total width of the corresponding quarkonium. 
The diagrams for two-gluon decays of quarkonium can be easily obtained from the diagrams for two 
photon decays of psendoscalar quarkonia (in Fig.5), with a simple replacement of photons by gluons, and 
hence, the two quark-photon vertices with the corresponding quark-gluon vertices. This would lead to the 
replacement: a —)■ IH]; in the expression for Fp in Eq.(58), entering into the two-photon decay 

width formula, with as being the QCD coupling constant. Here we have taken = -|-|-\/ 47 ra for c-quark, 
and Cq = — t-y/dvra for b-quark. The results are shown in Table 9. 
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BSE-GIA 

Expt. 

BSE[29] 

Pot. Model |2T] 

BSEjH] Pot. Model[l6| 


13.0699 

26.7± 3.0 

19.6 

32.44 

10.57 9.010 

^r]c{2S)^gg 

9.5340 

14.0± 7.0 

12.1 

24.64 

5.94 

^r]c{3S)^gg 

4.4123 



21.99 


^Vb{3S)^gg 

10.8646 


6.98 

13.72 

12.39 

^Vb(‘2S)^gg 

7.4766 


3.47 

6.73 

5.61 

^rib{3S)^gg 

6.2763 



5.58 

4.11 


Table 9: Two-gluon decay widths of ground state (IS) and excited state (2S) and (3S) pseudoscalar mesons, 
rjc and r/h in present calculation (BSE-CIA) along with their masses in other models and experimental data 
(all values are in units of Mev). 

6 Numerical Results and Discussions 

We have employed a 3D reduction of BSE (with a 4 x 4 representation for two-body {qq) BS amplitude) 
under Covariant Instantaneous Ansatz (CIA) for deriving the algebraic forms of the mass spectral equations 
whose analytic solutions (both eigen functions and eigen values), in Eq.(42-43), lead to mass spectra for 
ground and excited states of both pseudoscalar (r/c, and rjb) and vector (T/T, and T) quarkonia, in an 
approximate harmonic oscillator basis. The masses, and the eigen functions so obtained are used for 
calculating the leptonic decay constants, weak decay constants, two-photon decay widths, and two-gluon 
decay widths for ground and excited states of these pseudoscalar and vector quarkonia. 

The mass spectrum calculated in this BSE framework for (IS*,...,45') states of rjc and rjb, while for 
(IS*, 2S, ID, 3S, 2D,...,5D) states of J/'ijj and T are shown in Table 3 and 4 respectively. All numerical 
calculations have been done using Mathematica. We selected the best set of 6 input parameters, that gave 
good matching with data for masses of ground and excited states of rjc, rjb, J/T and T mesons. This input 
parameter set was found to be Co = 0.21, a;o=0.15 GeV., A=0.200 GeV, and Ao=0.01, along with the 
input quark masses rric = 1.490 GeV. and rrib = 5.070 GeV. The same set of parameters above was used to 
calculate the leptonic decay constants of ric,r]b, J/'ip, and T, as well as the two-photon and two-gluon decay 
widths of rjc, and rjb. However the experimental data on masses and decay constants/ decay widths of 
many of these states is not yet currently available. The results obtained for masses of ground and radially 
excited states of r]c, r]b, shown in Table 3, are in reasonable agreement with experiment. However, a wide 
range of variation in masses of various states in different models such as Lattice QGD model [8] , and QGD 
sum rule model [12] can be seen in Table 3. 

As regards the mass spectral predictions for vector mesons is concerned, we have listed the values for 
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Mj/^ for IS", 2S', 1D,...5S', and for Mt for IS", 2S', in Table 4. Many other states snch as SD, 

4Zi), and 5D for J/T are not yet experimentally available. The same holds trne for all the D states of 
T. For vector cc qnarkonia, the masses of gronnd and excited states are very close to central valnes of 
data. However the 15, and 25 states of bb are somewhat over estimated from central valnes of data. The 
disagreement with data increases as we go to the higher excited states of bb. This is dne to the fact that we 
did not incorporate the one-glnon-exchange (OGE) effects in the kernel, and nsed only the conhning part 
of interaction as in onr Eq.(31) (taking analogy from [191 ESI ESI SB] for heavy qnarkonia). Onr resnlts 
reflect the fact that the OGE term becomes more and more important as we go to very heavy bb qnarkonia, 
where the distance between the qnark and the anti-qnark will be tiny, and its contribntion to their mass 
spectra will be snbstantial. Farther, there is a degeneracy in the masses of 5 and D states with the same 
principal qnantnm nnmber N for J/T and T. The inclnsion of OGE terms in the potential will also lift 
np the degeneracy in these states and also bring masses of bb vector qnarkonia closer to experiment. 

However, we wish to mention that in this paper, onr main emphasis was to show that this problem 
of 4 X 4 BSE nnder heavy qnark approximation can indeed be handled analytically for both masses, as 
well as the wave fnnctions in an approximate harmonic oscillator basis. The analytical forms of wave 
fnnctions obtained as solntions of mass spectral eqnations (that are derived from 3D BSE), were then 
used to calculate the decay constants and decay widths for various processes involving these qnarkonia. 
We next intend to incorporate the OGE effects perturbatively into our mass spectral equations, using our 
wave functions so obtained in Eqs.(41) as unperturbed wave functions in this study. 

As regards our results for decay constants is concerned, it is seen that our calculated fp value of ground 
state is fr]c{is) = 0.4044GeV., and is within the error bars of its experimental value, 0.335 ±0.075Gev [34] . 
However a wide range of predictions of decay constant values for all states can be seen in different models 
in Table 6. The experimental data for many of these states is not yet available. However, it is observed 
that the decay constants keep decreasing as one goes from (15) to (35) states for both rjc and rji, mesons. 
This also signihes that the instability of these states increases with increase in radial quantum number N. 
This trend is similar to the trend observed in a recent potential model calculation [21] of decay constants 
of Tjc and rjh. 

Our results for leptonic decay constants of ground and excited states of J/V’ and T are listed in Table 
7. Our decay constants for J/'ij) and T for excited states are somewhat on the higher side in comparison 
to central values of data for these states. However, the decay constant fy values of various models again 
show a very wide range of variation as can be seen from Table 7. 

Then, with the same set of parameters, we calculate the two-photon radiative decay widths of the 
ground and hrst two radially excited states of r]c and rjb mesons. The two photon decay widths of our model 
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are listed in Table 8. Our result for two-photon decay width of ground state of rjc is r^^(lS') = 7.9178KeV., 
and is within the error bars of data, Tn^^is}{Exp.) = 7.2 ±1.2 KeV [35]. We have compared our results 
with those of other models, though data on many of these is not yet available. The two photon decay 
widths of other models also show a wide range of variation. 

We then calculate the two-gluon decay widths in our model for ground and radially excited states of 
r]c and rjf,, which are shown in Table 9. The two-gluon decay process accounts for a substantial portion 
of hadronic decay widths for states below cc or bb threshold. However, as pointed out in m, due to 
significant contributions from radiative corrections as well as from three-gluon decays, the two gluon mode 
does not give the complete picture. Our results on two-gluon widths are thus smaller than the hadronic 
widths of rjc and rji, states. However as can be seen from Table 9, the results of two-gluon decay widths in 
various models again show a wide range of variations. 

However, as mentioned earlier, our main emphasis in this paper was to show that this problem of 4 x 4 
BSE under heavy quark approximation can indeed be handled analytically for both masses, as well as the 
wave functions for / = 0, and I = 2 states. The validity of heavy quark approximation for quarkonia is due 
to the fact that not only the relative momentum between heavy quarks in the bound states is considered 
small, but also these quarks are treated as almost on mass shell [53|, which is justihed for calculation of 
low energy properties like the mass spectrum, and the decays of quarkonia (cc and bb systems). With 
the above approximation, that is quite justihable and well under control, in the context of heavy quark 
systems, we have been able to give analytical solutions of mass spectral equations Eqs. (35-36) of both 
pseudoscalar and vector quarkonia, giving us a much deeper insight into this problem. 

Analytical forms of the wave functions for both S and D states for n = 0,1, 2,... for cc, and bb systems 
thus obtained are given in Eq.(41). These wave functions were then used to calculate the decay constants 
and decay widths for various processes involving these quarkonia. We have also plotted the graphs of all 
the wave functions for the states, IS", ...,45" for rj^ and rji, mesons, and for the states IS", 2S, ID, 3S, 2D, 45", 
and 3D for J/T and T mesons, in Figs. 1-4. The over all features of all these plots show that the states 
nS and nD have n — 1 nodes in their wave functions. 

We are not aware of any other BSE framework, involving 4x4 BS amplitude, that treats the mass 
spectral problem involving heavy quarkonia analytically. To the best of our knowledge, all the other 4x4 
BSE approaches treat this problem numerically just after they obtain the coupled set of Salpeter equations 
(see [H]). We further wish to mention that the over all features of our plots of wave functions in Eqs. (41) 
derived in our framework are very similar to the corresponding plots of wave functions obtained by purely 
numerical methods in [18] , suggesting that our approach is not only in good agreement with the numerical 
approaches followed in other works, but also gives a deeper insight into the problem, by showing an explicit 


dependence of the mass spectrum on principal quantum number N as in Eqs. (42-43). As mentioned above, 
we next intend to incorporate the OGE effects into the mass spectral equations perturbatively. However, 
this we intend to do as further work. 

We also intend to extend this study to calculation of observables of heavy-light mesons such as D 
and B, which would involve incorporation of unequal mass kinematics in our framework, which is beyond 
the scope of the present paper, and we wish to do as further work. This study will also be extended to 
processes involving quark-triangle diagrams with two or more hadronic vertices such as to decays V —>■ P 7 , 
and V —>■ PP (with P and V being the pseudoscalar and vector quarkonia respectively), as further works. 
It is further expected that our present framework (with 4x4 representation of two-body BS amplitude) 
will be able to overcome the complexities in the amplitudes for these processes that appeared in our earlier 
framework |42j. 
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